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PREFACE 


This  Memorandum,  prepared  for  the  Advanced  Research  Projects  Agency, 
is  part  of  a  study  of  those  phenomena  which  affect  the  performance  of 
optical  or  infrared  reconnaissance  and  guidance  equipment*  The  objec¬ 
tive  of  these  studies  is  to  provide  sufficient  understanding  for  the 
system  analyst  to  compute  performance  estimates  under  various  opera¬ 
tional  conditions. 

A  quantitative  understanding  of  the  effect  of  atmospheric  turbu¬ 
lence  in  reducing  the  lateral  coherence  of  an  initially  coherent  wave- 
front  is  required  for  the  prediction  of  the  performance  of  various 
devices  employing  lasers  for  target  acquisition  or  guidance  in  tac¬ 
tical  missions.  Such  applications  are  characterized  by  near-horizontal 
propagation  paths  near  the  ground  on  the  order  of  one  to  tans  of  kilo¬ 
meters  in  length. 

These  results  should  be  of  use  to  those  interested  in  tactical 
applications  of  laser  range  finders,  laser  line  scanners,  and  the  vari¬ 
ous  guidance  systems  employing  an  illuminating  or  focused  beam. 
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SUMMARY 


It  is  shown  that  the  most  commonly  used  expression  for  the  mutual 
coherence  function  for  an  optical  wave  propagating  in  a  turbulent  at¬ 
mosphere,  based  on  an  unphysical  extrapolation  of  the  Kolmogorov  spec¬ 
trum,  is  incorrect.  In  modifying  the  spectrum,  it  is  shown  that  the 
corrections  to  the  coherence  function,  the  implied  resolution,  spread¬ 
ing  of  a  finite  beam,  and  the  signal-to-noise  ratio  using  heterodyne 
detection,  can  be  considerable.  Approximate  expressions  for  the  co¬ 
herence  function,  valid  over  three  distinct  propagation  distance  re¬ 
gimes,  are  derived.  In  general,  it  is  shown  that  the  field  has  a 
greater  transverse  coherence  length  than  that  predicted  by  other  authors, 
with  the  difference  being  more  pronounced  at  shorter  ranges.  This  im¬ 
plies  that  under  many  conditions  the  degradation  in  the  performance 
of  devices  which  depend  on  the  above  properties  is  less  than  previously 
reported . 
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I.  INTRODUCTION 


The  mutual  coherence  function  (defined  as  the  cross-correlation 
function  of  the  complex  field  in  a  direction  transverse  to  the  direc¬ 
tion  of  propagation)  is  the  quantity  which  describes  the  loss  of  co¬ 
herence  of  an  initially  coherent  wave  propagating  in  a  turbulent  medium. 
As  a  result »  the  mutual  coherence  function  is  important  for  a  number 
of  practical  applications.  It  determines  the  signal-to-noise  ratio  of 
an  optical  heterodyne  detector*  the  limiting  resolution  obtainable  along 
an  atmospheric  path*  and  the  mean  intensity  distribution  of  an  initially 

coherent  wave  emanating  from  a  finite  aperture.  This  study  demonstrates 

(1-3) 

that  the  most  commonly  used  expression  for  the  mutiial  coherence 

function  (given  in  Eq.  (5)  below)  is*  in  general*  incorrect.  This  ex- 

(2) 

pression  is  derived  from  the  assumption  that  the  Kolmogorov  spectrum 
of  index  of  refraction  fluctuations  can  be  extrapolated  to  arbitrarily 
small  wave  numbers*  K*  for  the  purpose  of  computing  the  coherence  func¬ 
tion.  The  apparent  justification  for  using  the  Kolmogorov  power  law 

for  K  ^  L  ^  (where  L  is  the  outer  scale  of  turbulence)  is  that*  al- 
o  o 

though  the  extrapolated  spectrum  diverges*  the  integrals  necessary  to 
compute  the  coherence  function  from  the  spectrum  remain  convergent. 
However*  the  sensitivity  of  the  result  to  the  divergent  spectrum  proves 
to  be  considerable. 

Section  II  argues  against  the  use  of  an  unbounded  spectrum  (e.g.* 
the  extrapolated  Kolmogorov  spectrum)  and  demonstrates  the  insensitivity 
of  the  coherence  function  to  any  boundea  spectrum.  Based  on  these  con¬ 
siderations*  a  modified  von  Karman  spectrum  (which  levels  off  for 

K  £  L  ^)  is  examined*  and  the  insensitivity  to  a  number  of  'Reasonable'* 
o 

spectra  in  computing  the  coherence  functions  is  demonstrated.  Compari¬ 
son  with  the  previously  accepted  result  predicts  an  improvement  in  the 
implied  resolution*  with  greater  improvement  over  shorter  paths. 

Section  III  shows  that  there  exist  three  distinct  propagation  dis¬ 
tance  regimes  for  which  approximate  expressions  for  the  coherence  func¬ 
tion  can  be  found.  These  formulas  with  the  respective  ranges  of  validity 
are  presented  in  the  table  on  p.  8  for  the  modified  von  Karman  spectrum 
of  Section  II. 
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In  Section  IV  we  calculate  the  Implications  with  regard  to  reso¬ 
lution  and  beam  spreading  of  our  expression  for  the  mutual  coherence 
function  or  modulation  transfer  function  (KTF)  for  the  case  of  homo¬ 
geneous  Isotropic  turbulence.  In  particular,  we  show  the  coherence 
length  (defined  as  the  transverse  separation  at  which  the  KTF  la  equal 
to  e  la  generally  greater  than  the  previously  accepted  value. 

Finally,  In  Section  V  ve  examine  the  Implications  for  coherent 
optical  detection.  Comparing  our  results  for  heterodyne  detection 
with  those  of  Fried, ve  predict  a  greater  long-term  average  algnal- 
to-nolse  ratio.  Tn  particular,  where  Fried  finds  a  maximum  useful  re¬ 
ceiver  diameter  for  all  ranges,  ve  find  that,  for  distances  which  are 
small  compared  with  the  mean  field  decay  length  (defined  In  Section 
II),  the  slgnal-to-nolse  ratio  Increases  Indefinitely  with  receiver 
size. 
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II.  THE  SPECTRUM 


The  analysis  Is  based  on  the  expression  for  the  mutual  coherence 

function  or  MTF  for  the  case  of  a  plane  wave  Incident  upon  a  homogeneous, 

(2) 

Isotropic  turbulent  medium: 


M(p,z) 


(1) 


* 


where  k  Is  the  optical  wave  number,  p  Is  the  transverse  separation  at 
propagation  distance  z,  and  Is  the  three-dimensional  spectral 

density  of  the  Index  of  refraction  fluctuations.  Equation  (1)  may  be 
written  In  the  form 


M(p,z) 


1 00 

J  (Kp)  ♦  (K)K  dK 
iQ  o  n 


>00 

,  ♦  (K)K  dK 

Iq  n' 


(2) 


where 


z 

c 


(3) 


(2) 

can  be  shown  to  be  the  propagation  distance  In  which  the  mean  field 

of  a  plane  (or  spherical)  wave  (u)  decays  to  e  from  Its  vacuum  value 

(the  angular  brackets  denote  an  ensemble  average).  From  Eq.  (2),  It 

”*2z/  z 

follows  that  as  p  -►  M(p,z)  e  This  expression  Is  in  ac¬ 

cord  with  the  physical  picture  of  the  light  arriving  at  the  two  points 
r^y  ^2  being  scattered  through  statistically  independent  media  when 
^  *  III  "  £2^  sufficiently  large  that  M(p,z)  *  U(ll)U*(r2)  -► 

<()(Xl)>  <^*(12)) . 


In  deriving  Eq.  (1)  It  has  been  assumed  that  <1>  Is  not  a  function 

n 

of  propagation  distance.  The  modification  of  Eq.  (1)  to  Include  an  ex¬ 
plicit  dependence  on  range  Is  z<t'^(K)  -*•  /q  <t'jj(K,z^)  dz^  .  For  spherical 
wave  propagation  z[l  -  J^(Kp)  ]  ->■  J*  dz'  [1  -  J^(Kpz'/z)]. 


(2) 

The  Kolmogorov  spectral  density  Is  most  comnonly  used  to  repre¬ 
sent  atmospheric  Index  of  refraction  fluctuations: 


^  (K)  -  0.033C^k"^^^^ 
n  '  n 


(A) 


It  Is  valid  within  the  Inertial  subrange  <<  K  «  ■fj't  where  e 

2vJC  and  L  e  2it£  are  the  Inner  and  outer  scales  of  turbulence,  re- 
o  o  o 

spectlvely,  and  Is  the  Index  structure  constant. 

In  order  to  compute  the  MTF  from  Eq.  (1),  It  Is  necessary  to  make 

certain  reasonable  assumptions  regarding  the  spectrum  outside  of  the 

(2)  (1  3) 

Inertial  subrange.  Tatarskl  and  others,  '  using  the  spectrum  of 
Eq.  (A)  for  all  K,  have  cosiputed  an  HTF  given  by 


M^(p,z)  -  exp 


1  «  p  «  L 

0  0 


(5) 


which  Is  used  for  all  z,  the  apparent  justification  being  that  the  Inte¬ 
gral  In  Eq.  (1)  converges.  However,  the  sensitivity  of  the  MTF  to  the 
extrapolation  of  Eq.  (A)  for  K  <  £  ^  Is  considerable.  This  can  be  seen 

®  _i 

by  carrying  out  the  Integration  In  Eq.  (1)  for  <  K  <  <»,  which  gives 
M'  (p,z)  -  exp  |-  kVzp^^^  [l  -  0.67(p/t^)^^^  +  0(p/li^)^]  | 


for  I  «  p  «  L 
o  o 


(6) 


From  Eqs.  (5)  and  (6) 


M'  -  M 


Evaluating  Eq.  (7)  at  that  value  of  p 
M  Is  equal  to  e  ^  yields 


(7) 


where 
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2  -15  -2/3 

For  example,  with  X  ■  10.6  u,  C  ■  3  x  10  cm  ,  L  ■  100  cm  (in 

n  o 

this  case  z  ^  1  km) ,  we  obtain  a  percentage  error  In  the  MIF  of  38  per- 
c 

cent  for  z  ■  1  km  and  23  percent  for  z  ■  10  km. 

This  example  reveals  the  sensitivity  of  the  KTF,  In  the  Inertial 
subrange,  to  the  physically  unreasonable  extrapolation  of  Eq.  (A)  for 
K  <  Although  the  Integral  In  Eq.  (1)  converges,  an  extrapolation 

would  lead  to  divergent  Integrals  for  both  the  energy  per  unit  volume 
of  the  fluctuations  and  the  distance  over  which  the  mean  field  decays 
(l.e.,  .  On  physical  grounds,  the  spectrum  must  begin  leveling  off 

for  K  corresponding  to  scales  large  compared  with  the  separations  over 
which  the  temperature  fluctuations  exhibit  appreciable  correlation, 
with  a  finite  upper  bound  as  K  0.  Further,  careful  examination  of 
Eqs.  (1)  or  (2)  reveals  the  Insensitivity  of  the  Integrals  to  the  spec¬ 
tral  density  for  any  spectrum  tdilch  remains  bounded  for  K  <  and 

hence,  for  the  purpose  of  computing  the  KTF,  we  suggest  the  use  of  any 

bounded  spectrum  In  this  range.  The  spectrum  falls  off  very  rapidly 

-1  (A) 

for  K  >  due  to  viscous  damping,  and  It  Is  customary  to  use  a 

gausslan  decay  In  this  region. 

An  example  of  a  spectrum  which  Is  convenient  for  computational 
purposes  Is  the  modified  von  Katman  spectrum 


♦n(K) 


0.033C„e  ° 

o 


(8) 


which  Implies  a  flat  spectrum  for  K  <  2.^^.  For  example,  substituting 
Eq.  (6)  Into  Eq.  (3)  yields,  for  «  L^, 


z 

c 


(9) 


Figure  1  Is  a 
wavelength  X. 


graph  of  the  mean  field  decay  length  z^  plotted  versus 

for  2:  ■  100  cm  and  typical  values  of  C^. 

o  n 
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III.  APPROXIMATE  FORMULAS  FOR  THE  MTF 


While  It  is  a  simple  numerical  calculation  to  compute  the  MTF  di¬ 
rectly  using  the  spectrum  of  Eq.  (8)  in  Eq.  (1),  It  Is  useful  to  have 
approximate  formulas  for  estimating  the  coherence  length  at  various 
ranges .  Based  on  Eq.  (8) ,  we  derive  approximate  expressions  for  the 
MTF  valid  over  particular  propagation  paths. 

First,  considering  the  Inertial  subrange,  substituting  Eq.  (8) 

Into  Eq.  (1),  assuming  p  »  1^,  and  expanding  to  lowest  order  In  (p/£^), 
we  obtain  (sae  the  appendix) 

Mj^(p,s)  -  exp  k^cjsp^^^j^l  -  0.80(p/i^)^^^]|  (10) 

Comparing  Eq.  (10)  with  Eq.  (6)  demonstrates  the  Insensitivity  of  the 
MTF  to  the  form  of  the  spectrum  for  K  <  provided  It  ramslns  bounded 
for  K  -»■  0.  At  a  given  range,  for  all  of  tha  transverse  separations  of 
Interest  to  lie  In  tha  Inertial  subrange  It  Is  necessary  that  M(£  ,s)  s*  1 
and  M(L^,z)  «  1.  This  Is  essentially  tha  condition 

z  «  X  «  z.  (11) 

c  1 

where  z^  Is  given  by  Eq.  (9),  and  z^,  obtained  by  replacing .2;^  *^0 

In  tha  formula  for  z^.  Is  the  dlstsnce  at  which  the  coherence  length 
of  the  field  Is  of  the  order  of  the  Inner  scale. 

For  ranges  greater  than  z^^,  all  of  the  p's  of  Interest  are  small 
compared  with  tha  Inner  scala.  The  Bassel  function  In  Eq.  (1)  can  then 
be  expanded  In  powers  of  p/C^  to  yield 

M(p,z)  a-  exp  1.72k^C^r^^^zp^^  z  »  z^  (12) 

Finally,  for  z  «  z^,  recalling  M(»,z)  -  »  i,  the  MTF  Is 

essentially  unity  for  all  values  of  p,  l.a.. 


z  «  z 

c 


M(p,z)  -  1 


(13) 


-8- 


In  Fig.  2  we  compare  the  KTF  computed  from  Eq.  (5)  (curve  labeled  M  ) 

o 

with  that  obtained  from  Eq.  (1)  using  the  spectrum  given  by  Eq.  (8) 
(curve  labeled  M)  and  Eq.  (10)  (curve  labeled  Mj^) ,  taking  X  •  10.6  w, 
z/Zj,  ■  10,  ■  1  ®ni»  and  ■  2it  m.  In  genei.‘al,  the  present  analysis 

Indicates  a  "more  coherent"  MTF. 

The  approximate  formulas  for  the  plane  wave  MTF  with  the  respective 
ranges  of  validity  are  sunmarlzed  In  the  table. 


APPROXIMATE  EXPRESSIONS  FOR  THE  MODULATION  TRANSFER  FUNCTION 


Range 

MTF 

z  «  z 

c 

1 

exp  kVzp^^^[l  -  0.80(p/L^)^^^1j-  - 

Z  Z  Z  a 

C  1 

1 

exp  {-  3.72(z/z^)(p/«.^)^^^[l  -  0.80(p/«.J'/^^]} 

z  »  z^ 

f  ,  „.2„2  -1/3  2\ 

exp  •<-  1.72k  zp  j  ■ 

exp  2.4(z/zj^)(p/£^)^| 

NOTE:  The  quantities  z  (the  distance  where  the  average  field  Is 
-1  ^ 

down  by  e  )  and  z^  (the  distance  where  the  coherence  length  of  the 
field  Is  of  the  order  of  the  Inner  scale  of  turbulence)  are  given  by 
(0.39k^cV^^)"^  and  (o. 39k^cV^^)"^,  respectively. 


4 


uojioun^  j8jSU0J4  uoj40|npo^ 


Compari$on  of  MTF  as  a  function  of  p 
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IV.  RESOLUTION  AND  BEAM  SPREADING 


There  are  two  Important  consequences  of  our  "more  coherent"  MTF. 

The  mean  visibility  Is  better  than  predicted  by  Frled^^^  for  given 

turbulence  parameters  along  a  uniform  atmospheric  path,  and  the  average 

amount  by  which  a  finite  beam  spreads  Is  smaller. 

Defining  as  the  transverse  separation  at  which  the  atmospheric 

MTF  Is  reduced  by  e  the  minimum  resolvable  length  at  a  distance  z 

from  an  observer  Is  well  known  to  be  ~  z/kp  .  It  can  be  shown  that 

o 

when  p  Is  small  compared  to  the  size  of  the  transmitting  aperture, 

®  f  61 

the  angular  spread  of  a  finite  beam  due  to  the  atmosphere  Is  ~  1/kp^.^  * 

It  follows  from  the  results  of  the  previous  section  that,  for  dls- 

5  3 

tances  small  compared  with  z.(  >  10  z  for  L  /l  >  10  ) ,  the  plane 

1  c  o  o 

wave  MTF  does  not  depend  on  the  Inner  scale,  and  can  be  written  as 


M(p,z) 


f 


Jo  [U  + 


Jq(u)u  du 


■  F(p/^  ,z/z  )  z  «  z. 

O  C  1 


(14) 


The  modification  of  Eq.  (14)  for  spherical  waves  Is  given  by 


M  (p,z)  -  exp 


(  du  uJ^(u)  ^  ^ 

^0  Jo  [u  +  (solZ  ) 


(p/i  )^^^  \  du  uJ  (u) 
o  1«  o 


S  F 


z/z  ) 
c 


(15) 


In  Figs.  3a  and  3b  we  plot  p  /-L  versus  z/z  for  plane  and  spherical 

o  o  c 

waves,  which  are  obtained  by  Inverting  the  equations  F(p  ,  z/z  )  >•  e 

o  o  c 

and  T  (p  /Z  ,  z/z  )  ■  e  ,  respectively.  On  the  same  graph  In  each  fig- 
s  o  o  c 

ure,  we  compare  our  results  with  that  obtained  using  the  MTF  of  Eq.  (5) 
(dashed  curve)  and  the  corresponding  expression  for  spherical  waves 
(obtained  by  multiplying  the  exponent  In  Eq.  (5)  by  3/8).  For  the  pur¬ 
poses  of  this  computation,  we  have  written  Eq.  (5)  In  the  form 
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M^(p,z)  -  exp 

The  percentage  error  in  each  case,  defined  as  I (p'  -  p  )/p  I  x  100 

0  0  0 

where  p*  satisfies  M  (p^ ,z)  ■  e  ,  is  indicated  on  the  scale  at  the 

O  0  0 

right. 

As  one  approaches  the  plane  or  spherical  wave  source,  the  present 
analysis  indicates  that  the  transverse  coherence  of  the  field  Increases 
at  a  greater  rate  than  previously  predicted.  The  previous  results  are 

A 

in  error  by  ~  10  percent  at  10  z  and  100  percent  at  0.5z  ,  while  for 

distances  <  0.5z  ,  the  "new"  MTF  is  never  down  to  e  .  For  z  >  z. 

5  3  ^ 

(■  10  z^  for  ■  10  ) ,  the  use  of  again  gives  a  poor  approxi¬ 

mation  to  the  "new"  MTF,  which  in  this  region  is  given  by  Eq.  (12). 


-13- 


V.  SIGNAL-TO-NOISE  RATIO  IN  HETERODYNE  DETECTION 

As  an  additional  illustration  of  the  Implications  of  our  "more  co¬ 
herent"  MTF,  we  use  it  to  compute  the  signal-to-nolse  ratio  in  a  coherent 
detection  system  and  compare  the  results  with  those  computed  by  Fried. 

The  slgnal-to-noise  ratio  is  given  by^^^ 

-  4(n/e)(A^)D^  £  xK^(x)M(Dx,z)  dx  (16) 

where  is  the  signal  amplitude,  n/e  is  the  quantum  efficiency  measured 
in  electrons  per  unit  energy,  D  is  the  diameter  of  the  collecting  aper¬ 
ture,  and 

Ko(x)  ■  J  [co«"^(x)  -  x(l  -  (17) 

Then  Eq.  (16)  can  be  %n:ltten  In  the  form 

A(n/e)  (18) 

for  0  <  z  «  Zj^,  where 

'!'(A  t/t)  xK^(x>F(5x,  t/t)  dx  (19) 

c  Jo  °  c 

The  function  V  contains  the  dependence  of  slgnal-to-noise  ratio  on  the 
collector  diameter  in  the  presence  of  atmospheric  distortion  and  is 
plotted  as  a  function  of  ^  in  Figs.  4a  and  4b  for  various  values  of 

z/z  .  The  reduced  slgnal-to-noise  ratio  Y,  as  computed  from  Eq.  (19) 

®  (1) 

(solid  curve),  is  compared  with  the  results  of  Fried, ^  *  who  used  the 

MTF  given  by  Eq.  (5)  (dashed  curve).  In  general,  the  present  analysis 

predicts  a  somewhat  larger  signal-to-nolse  ratio.  The  difference  is 

more  pronounced  for  distances  small  compared  with  z^,  where  (s)/N  as 

derived  here  Increases  indefinitely  with  aperture  size.  This  result 

is  in  contrast  with  that  of  Fried,  where  an  effective  limiting  diameter 
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Is  derived  beyond  which  Increasing  the  diameter  results  in  very  little 
Improvement  in  {s)/N  for  all  ranges.  For  ranges  very  large  compared 
with  (but  still  small  compared  with  ,  the  MTF  of  Eq.  (1)  approaches 

that  of  Eq.  (5),  and  the  difference  between  the  two  slgnal-to-nolse 
ratios  tends  to  zero.  This  trend  Is  Illustrated  In  Fig.  5,  where  the 
reduced  slgnal-to-nolse  ratio  Is  plotted  as  a  function  of  z/z^  for  ■  1. 
The  solid  curve  Is  computed  from  thfl  MTF  given  by  Eq.  (14)  while  the 
dashed  curve  Is  computed  from  the  MTF  of  Eq.  (5).  The  percentage  error, 
defined  as  the  difference  between  {s)/N  using  Eq.  (1)  and  Eq.  (5), 


Fig.  5 — Reduced  signal -to-nolse  ratio  at  10.6  ^  in  a  heterodyne  detection 
system  as  a  function  of  for  diameter  of  the  collector  equal  to 
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divlded  by  (s)/N  using  Eq.  (1),  multiplied  by  one  hundred,  Is  also 
plotted  In  Fig.  5  (referring  to  the  right-hand  ordinate).  We  see  that 
a  maximum  percentage  error  of  53  percent  Is  obtained  for  k.  In 

general,  the  range  corresponding  to  maximum  percentage  error  will  of 
course  vary  with  the  ratio 
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Appcndlx 

EVALUATIOW  OF  AN  INTEGRAL 


Consider  the  Integrel 


KPiO 


i0.033cj  exp  [-  (K*)^J  ) 


K  dK 


(A-1) 


obtained  by  eubetltutlng  the  epectrue  given  by  Eq.  (8)  Into'Eq.  (1). 
Subetltutlng  x  ■  Kp  In  Eq.  (1)  ve  obteln 


I(p,i)  ■  4wVi(0.033)C^p*^^ 

n 


-  (1  -  J  (x)]  exp  I-  (xE  /p)^l 


7777^ 


X  dx  (A-2) 


where 


a 


2 


(A-3) 


For  thoee  veluee  of  p  in  the  inert  lei  eubrenge  the  exponential  in 
Eq.  (A-2)  can  be  neglected,  eince  the  xain  contribution  to  the  inte¬ 
gral  coxae  for  x  ^  p/t  «  p/i  .  Then  euppreeeing  the  factor 
4x^k^s(0.033)C:p^'  ,  we  obtain 


J-  [1  -  J  (x))x  dx 

2  2  11/6 
0  (x^  + 


1  «-5/3 

5  “ 


■t 


r(ll/6)a‘ 


(A-A) 


where  Is  the  nodlfied  Beeeel  function  of  the  second  kind  of  order  v, 
and  r(x)  Is  the  ganna  function.  In  obtaining  the  second  tern  on  the 
rlgbt-hend  side  of  Eq.  (A-4) ,  we  heve  used  Eq.  (2),  p.  434  of  Ref.  7. 
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The  modified  Bessel  function  of  the  second  kind,  K^i  is  related  to 
the  modified  Bessel  function  of  the  first  klnd|  through  the  relation 


K 


V 


2  sin  TTV 


(A-5) 


Using  the  power  series  expansion 


I  (a) 


0^^ 

■  (f)''  S 


f  2  y  y  V  k 

Ca  /^) 


2/  +  k  + 


1) 


(A-6) 


in  Eq.  (A-5)  and  substituting  In  Eq.  (A-4)  with  v  ■  5/6  yields 


K5/6(a)  ^  TT 


(a/2)'^^*^  (o/2)^^^  (o/2)^^^  (a/2)^^^^ 

r(i/6)  r(776)  r(ii/6)  ’  r(i7/6) 


(A-7) 


Hence,  from  Eqs.  (A-4)  and  (A-7)  we  obtain 


S"  [1  -  Jp(x)] 

0^7777^’“”' 


1.12  -  0.90o 


1/3 


+  0  •  1 3o(  + 


(A-8) 


Replacing  the  factor  previously  withheld  yields  the  expression  given 
by  Eq .  (10) . 
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